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Zero-index-metamaterials (ZIM) have drawn much attention due to their intriguing properties and novel 
applications. Particularly, in a parallel plated ZIM waveguide system with defects, total reflection or 
transmission of wave can be achieved by adjusting the properties of defects. This effect has been explored 
extensively in different types of ZIM (e.g., epsUon-near-zero metamaterials, matched impedance ZIM, or 
anisotropic ZIM). Almost aU previous literatures showed that only monopole modes are excited inside the 
defects if they are in circular cylinder shapes. However, the underlying physics for excited modes inside 
defects is wrongly ignored. In this work, we uncover that additional modes could be excited by theoretical 
analysis, which is important as it wiU correct the current common perception. For the case of matched 
impedance zero-index metamaterials (MIZIM), the additional dipole modes can be excited inside the defects 
when total transmission occurs. Moreover, we also observe the same results in Dirac-cone-like photonic 
crystals which have been demonstrated theoretically and experimentally to function as MIZIM. For another 
case of epsUon-near-zero metamaterials (ENZ), we find that additional higher order modes (e.g., tri-pole) 
can be excited inside the defects when total transmission happens. Numerical simulations are performed to 
verify our finding regarding the additional modes. 

In the beginning of the new century, the first metamaterial was achieved to realize the function of negative 
refractive index'. After that, the research on metamaterials''"' has made great progress. Recently, attention to 
zero-index-metamaterials (ZIM)" '"'"'^ '''' has been extensive. For instance, matched impedance zero-index 
metamaterials (MIZIMs), epsUon-near-zero metamaterials (ENZ), anisotropic ENZ. By utilizing ZIM, some 
applications and devices with novel functionalities can be realized, such as squeezing wave energy'" tailoring 
wave front'""", realizing total transmission and reflection in ZIM"* waveguide bending'", enhancing radiation 
from an embedded source'''' controlling energy flux''', etc. Several years ago, by putting perfect electric con- 
ductor (PEC) or perfect magnetic conductor (PMC) defects in ZIM in a waveguide structure, Hao et aZ."* 
confirmed that incident electromagnetic wave can undergo total reflection or transmission. Later, Nguyen 
et aV^ found similar effects by introducing dielectric defects into MIZIM. However, due to an insufficient 
expression of the magnetic field, some interesting phenomena and physics are thereby missing. For example, 
they claimed that for total transmission, only monopole modes exist in the dielectric defects. A lot of publica- 
tions'"' "'-'^'*"' (including one from the authors^") later followed this erroneous step, albeit some other intriguing 
properties found. In this letter, we will give a more comprehensive analysis and show that additional higher modes 
are excited together with the monopole modes. Our paper corrects some common misunderstanding and shows 
more colorful physics for ZIM systems. 



Results 

Now let us start from the schematic plot of a two dimensional (2D) waveguide structure in Fig. 1. Region 0 and 
region 3 are free space. Region 1 is ZIM with the effective permittivity and permeability tii and ^i. Region 2 
consists of N cylindrical defects embedded in region 1. The effective permittivity and permeability of the j-th 
cylinder are e,2j and ^2j> respectively. Without loss of the generality, we suppose that a transverse magnetic (TM) 
wave (its magnetic field H is along z direction) is incident from the left port of the waveguide. The outer 
boundaries of the waveguide are set as PECs. If the incident wave is a transverse electric (TE) wave (with the 
electric field E polarized along z direction), the outer boundaries of the waveguide should be changed into PMCs 
for similar results. 

For simplicity, we assume the incident magnetic field Hint = Z-Hoze'''^°^^™'', where is the wave vector in free 
space with kg = cu/c, m is the angular frequency, c is the velocity of light in free space, Hqz is the amplitude of the 
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Figure 1 | The schematic description of the 2D waveguide structure. 

Region 0 and 3 is air. Region 1 is ZIM. Region 2 is the cylindrical defects. 
The parallel red lines are PEC boundaries of the waveguide. A TM wave is 
incident along x-direction from the left port of the waveguide. 

incident magnetic field. In the following sections, we will omit the 
time harmonic factor e^'"". The electromagnetic (EM) wave in each 
region satisfies the Maxwell's equations: 



V = . 



(1) 



where the integer m indicates each region and £„ is the relative 
permittivity of each region. The magnetic field in region 0 is a sum- 
mation of the incident wave and the reflected wave, and is written as. 



Ho = zHoz(e*°"' + Ke-*»^) 



the electric field is. 



(2) 
(3) 



where 5R is the reflection coefficient. 

Likewise, we can obtain the magnetic field and the electric field in 
region 3 as. 



(OEq 



(4) 



(5) 



where Q? is the transmission coefficient. In region 1, as almost 

equals to zero, V x Hi must be zero in order to guarantee a finite 
El. Consequently, the magnetic field in region 1 denoted as Hi should 
be a constant. By applying the boundary conditions at the interfaces 
of X = 0 and x = a, we have. 



+ l)Ho, = Hi and Hi = QHo, 



(6) 



Therefore, + 1 = 9. 

In region 2, the magnetic field inside each cylindrical defect fol- 
lows the Helmholtz equation, 

V'H2+e2,H2jklH2 = 0. (7) 



The solution can be written as a summation of infinite number of 
Bessel functions with angular terms. Therefore, the magnetic field in 
region 2 should be written as. 



_ N GO 

«2=zE E tj„J„{k2jrj)e'"' 



(8) 



where tj„ are the coefficients to be determined for the «-th order 
Bessel functions. 



By applying Dirichlet boundary conditions at the surface of each 
defect, the magnetic field in region 2 becomes. 



;=1 V-" 



+ aj„J„{k2jrj) cos{n6j) + Pj,J„{k2jrj) sin(«e;) ,(9) 



where }„{x) is the n-th order Bessel function with JJJcjjRj) = 0. aLj„ 
and ;6j„ are coefficients of the excited higher order modes, 
= ko^Jsajlhj is wave vector in the j-th cylindrical defect, Rj is 
the radius of the j-th cylinder, r, is relative radial coordinate in thej-th 
cylinder, Oj is relative angular coordinate in the j-th cylinder, as 
mentioned in Ref 19. We note that the original expression of the 
magnetic field in each defect is wrongly assumed in Ref 19 from the 
two missing terms that typify the additional modes in the defects. 
However, if JJJijjRj) 0, the two additional terms should not be 
included so that at each circular boundary r, = Rp the magnetic field 
takes the constant value Hi. 

With Eq. (9), the electric field inside each defect could be obtained 
as follows by recalling Eq. (1), 

^-^^^^ - =<;■„/„' (/C2,0) cos(«0,)-/}j.„7;(/:2,r,) sin(n0,) ) ,p0,,(lO) 



where Oj is the azimuthal unit vector for the j-th cylindrical defect. 
By using the Maxwell-Faraday equation. 



oEd\ = 



^ dt 



(11) 



and after some meticulous calculations, we could obtain the trans- 
mission coefficient as^°. 



1 



ikoiii{S-Sj,) 
2h 



N 

h 



(12) 



where S = a X h is the entire area of region 1 and region 2, 

N 

Sde = tiRJ is the total area of region 2, which consists of N 

i=i 
cylinders. 



For the system of MIZIM 

For matched impedance zero-index materials (MIZIM), fii '- 
(12) changes into the following formula. 



RjhikyRj) 



/o(M;) 



; 0, Eq. 



(13) 



which has been derived in Ref 19 (the Eq. (8) therein). From the Eq. 
(13), we see that to achieve total transmission (9— > 1), ]i{k.2jRj) must 
be equal to zero. Therefore, we should select « = 1 in Eq. (9), the 
magnetic field in region 2 is then written as. 



-zHiY. 



/o(M;) 



^ajih{k2jrj)cosej + PjJi{k2jrj)smej . (14) 



It is not only a function of but also a function of Bj, aji and (iji are 
coefficients to be determined. Eq. (14) is the unique solution for total 
transmission with JiikyRj) = 0. According to Eq. (14), we note that the 
magnetic field inside each defect consists of not only monopole modes, 
but also additional dipole modes. 

To solve the coefficients of dipole modes, we suggest an approx- 
imate system, i.e., a dielectric cylinder embedded in MIZIM as the 
background. As the waveguide system supports TMq mode, which 
resembles a plane wave, we suppose the incident wave is a plane wave 
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along x-direction. The magnetic field in MIZIM could be expressed 



H,^ = zHi J2 [i"J„{kir)+g„H,ikir)]e'"", (15) 

n — — 

where H„(x) is the n-th order Hankel function of the first kind, Hi 
takes the same value of the constant magnetic field in the MIZIM 
area, g„ is the scattering coefficients, ki = kgy^EiHi the wave vector 
in MIZIM. When total transmission occurs, the scattering coeffi- 
cients should be very tinny, that is g„ ~ 0. 

The magnetic field in the dielectric cylinder could be written as, 



H2=z ^ tij„(k2ir)e 



At the boundary of the cylinder r = R^, we should have. 
Hi f2 i"Jn{kir)e'"'= £ h„}„(k2ir)e'"" . 

n = — Oj n = —orj 

We can easily obtain the coefficient of each mode, 

UkiRi) 



(16) 



(17) 



tin — i"Hi 



Jn{k2lRl 



(18) 



By combining Eq. (16) and (18) with Eq. (9), we could obtain the 
relationship between ti„ and follows, 



ai„ = ti„/Hi=t ^ ^ , , n = 0, 



ai„=2ti„/Hi=i' 



In{k2lRl) 

„ 2J„{kiRi) 

Jn{k2lRl) ' 



«>1. 



(19) 
(20) 



For MIZIM, as Eq. (14) only has two terms (n = 0 and n = 1). We 
could therefore obtain that. 



— ^10 — 



an -■ 



JojkiRi) _ 1 
/o(fc2ii^i) /o(*:2ii^i)' 

. 2h(kiRi) 
'hikiiRiY 



(21) 



22 



ao is the coefficient of the monopole mode, and an is the coefficient 
of the dipole mode. Due to the symmetry in ^-direction, the degen- 
erate state of dipole mode inside the defect {J„{k2jrj)sm{nOj) in Eq. 
(9)) could not be excited out, the coefficient of the degenerate state 
should be zero. 

We plot the relationship between the parameter |o!ii|/|o:o| and 
different frequencies for three types of MIZIM (with Ei = ^tii = 
10"', 10-" and 10"=) from Eq. (21) and (22), as shown in Fig. 2. 
We also plot the related numerical result from the waveguide struc- 
ture (the simulation results in the following sections are all from 
COMSOL). We find that dipole mode is much more dominative than 
monopole mode at the resonance frequency that 7i(fc2i^i) ~ 0. 
While for the frequencies slighdy away from the resonance fre- 
quency, dipole mode becomes very weak or diminishes. When £i 
= fii tends to zero gradually, the coefficient of dipole mode will 
decrease accordingly, and the resonance peak becomes narrower. 
The resonance frequencies from the numerical results have tiny shifts 
due to the effect of PEC boundaries of the waveguide. When £i = 
tends to zero gradually, they will get closer to the analytic resonance 
frequency. 



We will verify the above findings from the numerical simulations. 
For simplicity, we assume that there is only one cylindrical defect in 
the region of ZIM. The radius of the cylindrical defect is 0.2 m. Its 
dielectric constant is £21 = 4, and its relative permeability is ^121 = 1. 
In order to achieve total transmission, the term 7i(fc2i^i) should be 
equal to zero. As a result, the working frequency is 0.45737 GHz. 
However, due to the effect of PEC boundaries of the waveguide, this 
frequency is slightly away from the real resonant frequency, leading 
to diminishing dipole mode. Therefore, in the simulations, we select 
the frequency of 0.457 GHz so as to obtain more clear dipole mode. 
In the region of ZIM, we set a = h = 0.8 m. The magnetic field 
distributions for the above waveguide system are shown in Fig. 3 
for ZIM with different permittivities. In Fig. 3(a), we set £1 = fii = 
0.001 and Hi = 1, and it seems that only monopole mode is excited 
from the field pattern. However, the field pattern just takes the real 
part of the magnetic field in COMSOL. In fact, if we read from the 
imaginary part, we could find that dipole mode exists. The magnetic 
field inside the defect is a summation of monopole mode and dipole 
mode from Eq. (14). The coefficients of monopole mode and dipole 
mode have a n/2 phase difference from Eq. (21) and Eq. (22). For 
instance, if we change Hi = 1 into Hi = i, we could find that dipole 
mode appears from the field pattern in Fig. 3(b). If we set Si = jii = 
0.0001, the dipole mode will become weaker, as shown in Fig. 3(c). 
For Ki = fii = 0.00001, dipole mode almost disappears, see in 
Fig. 3(d). We can also see this from Fig. 2. Numerically, we find that 
the coefficients of dipole modes are an = 12.37i, an = — 3.42i and 
an = —0.245! for the cases of Fig. 3 (b). Fig. 3 (c) and Fig. 3 (d), 
respectively. For all the above cases, the coefficients of the monopole 
modes are l/}o{k2iRi) = —2.483 and the coefficients of the degen- 
erate state of dipole modes are zero {pn = 0). For Ei = /.ii = 0.001, 
dipole mode is more dominative than monopole mode. While for £1 
= /ii = 0.0001, they are comparable to each other. For Ei = fii = 
0.00001, the monopole mode is more dominative than dipole mode. 
It seems that the dipole mode is diminishing when £1 = jii tends to 
zero gradually. It is not easy for us to choose the required resonant 
frequency as the resonance goes extremely narrow. For a particular 
value of near zero permittivity, we can in principle find a frequency 
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Figure 2\ agis the coefficient of monopole mode, an is the coefficient of 
dipole mode. Here we plot the ratio of the amplitude ofan and in log 
scale. The dash curves are numerical results, while the solid curves are 
analytical results. The black, red and blue curves are for the cases of £1 = jii 
= 10"-\ 10"'' and 10"^, respectively. The radius of the cylindrical defect is 
set to be 0.2 m. Its dielectric constant is £21 = 4, and its relative 
permeability is /Xji = 1- In the simulations, we set a = ?! = 0.8 m for the 
ZIM region. 
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Figure 3 | ThemagneticfielddistributionofMIZIM waveguide system with a dielectric defect, for (a) £i = = 0.001andHi = 1 (b) 8i = = 0.001and 
Hi = i; (c) £i = /ii = 0.0001 and Hj = i; (d) £i = /ij = 0.00001 and = i. The working frequency is 0.457 GHz. 



near the resonant one where dipole mode is much more dominative 
than monopole mode. 

To demonstrate the hybridization of monopole mode and dipole 
mode more clearly, we plot the real part of magnetic field distribution 
inside the defect from x = — 0.2 m to x = 0.2 m (2A y = 0) with 
different values of Hi for £1 = ^1 = 0.001, as shown in Fig. 4. For Hi 
= 1, only the information of monopole mode is observed from the 
black curve (it is an even function of x). While for Hi = i, the 
information of dipole mode could be observed from the red curve 
(it is an odd function of x). For Hi = 0.707 + 0.707i, both informa- 
tion of monopole mode and dipole mode could be observed from the 
blue curve. To make it more straightforward, we also plot the ampli- 
tude of the magnetic field (see the green curve), which is independent 
of Hi and has two symmetric peaks at the positions of x = 0.1 m and 
x= —0.1 m because of the existence of dipole mode. If there is only 
dipole mode inside the defect, the amplitude should be zero at the 
position of X = 0. However, the amplitude there is a value of about 
2.5, which is equal to the amplitude of the monopole mode at the 
position of X = 0 (see also the black curve). Therefore, both mono- 
pole mode and dipole mode exist inside the defect. 

As it is known, Dirac-cone-like photonic crystals"" can be regarded 
as MIZIM near Dirac point frequency. It should be possible to pro- 
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Figure 4 | The magnetic field distribution from x= —0.2 m to x = 0.2 m 

(at y = 0) in the defect. The black, red and blue curves are magnetic field 
distribution for the case of Hi = 1, Hi = i and Hi = 0.707 + 0.707;, 
respectively. The green curve is the amplitude of the magnetic field. 



duce the above similar effect if we replace MIZIM with such photonic 
crystals. It is noticed that the incident wave is now a transverse 
electric (TE) wave, and the outer boundaries of the waveguide struc- 
ture are PMCs. We will show that if a cylindrical defect is introduced 
in MIZIM or Dirac-cone-like photonic crystals, at the condition of 
Ji(k2iRi) = 0 when total transmission occurs, both dipole mode and 
monopole mode exist inside the defect in the waveguide system 
(^21 =fco \/£2iA'2i is wave vector of light in the defect, £21 and 
/.(21 are the permittivity and permeability of the defect respectively, 
Ri is the radius of the defect). Following Ref. 16, the Dirac-cone-like 
photonic crystals consist of cylindrical alumina rods arranged in a 
square lattice. The radii of the rods are 3.75 mm with a dielectric 
constant 8.8. The lattice constant is 17 mm. The Dirac point fre- 
quency/is about 10.3 GHz. We set a = /j = 0.187 m for the region 
of photonic crystals in the waveguide and insert a cylindrical defect 
with a radius of Ri = 0.0284 m in the center of the region (the radius 
should be large enough to visualize the above effect). The dielectric 
constant of the defect is £21 = 0.3905 and its relative permeability is 
/.(21 = 1, to satisfy Ji{k2iRi) = 0 at the Dirac point frequency. In Fig. 5 
(a), we plot the electric field for the system of Dirac-cone-like photo- 
nic crystals and choose a suitable phase of the incident plane wave (£1 
= /) so that only dipole mode is demonstrated in the cylindrical 
defect. Likewise, we choose another phase of the incident plane wave 
(£1 = 1), and plot the electric fieldin Fig. 5 (b), where only monopole 
mode is shown in the defect. For comparison, we replace the Dirac- 
cone-like photonic crystals with MIZIM, and plot corresponding 
electric field in Fig. 5 (c) and (d). Fig. 5 (c) shows a consistent dipole 
mode with Fig. 5(a), while Fig. 5(d) gives out monopole mode like 
that in Fig. 5(b). Therefore we cannot neglect the existence of dipole 
mode like Ref. 19. Sometime, it is more dominative than the mono- 
pole term near the resonance frequency, even for such a realistic 
photonic crystal system. 

For the system of ENZ 

After discussing the MIZIM case, we come to the ENZ case. Let us 
return to Eq. (12). In order to obtain total transmission, the following 
term must be zero. 



2h 



RjhikyR,) 



h{kyRj) 



(23) 



For simplicity, we suppose that there is only one cylindrical defect 
in the ENZ area. After some calculations, we shall have. 
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Figure 5 | The electric field distribution of Dirac-cone-Iike piiotonic crystals witli a defect for (a) £i = i and (b) £i = 1, and the electric field distribution 
of MIZIM waveguide system with a defect for (c) £i = i and (d) £i = 1. Here £i is the initial phase of the incident TE plane wave in the ZIM area. 
The radius of the defect is 0.0284 m, and its dielectric constant is 0.3905 and its relative permeability is /I21 =1. We set a = /i = 0.187 m for the ZIM 
area. The outer boundaries of the waveguides are PMCs. In addition, we set £1 = /Ji = 0.001 for the MIZIM in (c) and (d). The working frequency is 
10.3 GHz. 



Jl{k2lRl) k2lHi{S-Sde) 



Jo{k2iRi) 



2nRi 



(24) 



If the «-th order Bessel function /„{fc2i^i) is zero, the magnetic 
field in the defect could be written as. 



Joihiri) 
7o(fc2iRi) 



-^inJn{k2iri)cos{ne) + fi^J„{k2iri)sm(n9) (25 



Therefore, there is not only monopole mode excited inside the 
defect, but also some other additional higher order mode emerging 
as well, if both the conditions of Eq. (24) and J„{k2iRi) = 0 are 
satisfied. For example, we choose Jj,{k2iRi) = 0 and tune the config- 
uration and material parameters to satisfy Eq. (24). Following similar 
calculations to Eq. (19) and (20), we could get the coefficient of tri- 
pole mode 0!i3 = P2}j,{kiRi)/J^{k2iR]) and the coefficient of mono- 
pole mode 0!o = l//o(A:2i-Ri). Likewise, we find the relationship 
between the parameter |o(i3|/|ao| and different frequencies for three 
types of ENZ (with = 10"\ lO"'* and 10"=') both theoretically and 
numerically, as shown in Fig. 6. From the analytical results (solid 
curves), we find that tri-pole mode is much more dominative than 
monopole mode at the resonant frequency where 73(^21^1) = 0. For 
other frequencies slightly deviating from the resonant one, tri-pole 
mode is disappearing. When £] tends to zero gradually, the coef- 
ficient of tri-pole mode will decrease accordingly and the resonance 
peak will become narrower. However, the resonance frequencies 
from numerical results have a tiny shift because of the effect of outer 
PEC boundaries of the waveguide. When g] tends to zero gradually, 
the resonance frequency wiU approach the analytic resonance 
frequency. 

The finding will be confirmed again from numerical simulations. 
Suppose that there is one cylindrical defect inside the ENZ area. The 
radius of the defect is 0.01 m, its dielectric constant is g2i = 16 and its 
relative permeability is /(2i = 1. In order to make 73(^21^1) = 0, where 
fc2ii^i is the first root of the third order of Bessel function, the working 
frequency should be about 7.615 GHz. In addition, to satisfy Eq. (24), 



we set the effective permittivity and permeability of the ENZ as gj = 
0.001 and jii = 0.6 respectively and seta = h = 0.021 m for the ENZ 
area. The magnetic field distribution for the above system is shown in 
Fig. 7. The tri-pole mode is demonstrated inside the defect when Hi 
= i, as shown in Fig. 7(a). Likewise, we can also observe the mono- 
pole mode by changing Hi into 1, as shown in Fig. 7(b). We numer- 
ically find that the coefficient of tri-pole mode is c/.^ = — 1.67/, /^i3 = 
0 and the coefficient of the monopole is l/}o{k2iRi) = 4.17 (the tri- 
pole mode here is obvious but not dominative, to get a dominative 
tri-pole mode, the working frequency should be shift to about 



CO 



t/5 

o 

H-1 



-6 



io;V 

. - ' ' ' (V - ~ 


Num. 

Analy. 




















1 . 1 , 1 . 1 . 1 . 



7.6050 7.6075 7.6100 7.6125 7.6150 7.6175 7.6200 
Frequency(GHz) 

Figure 6 | otq is the coefficient of monopole mode, is the coefficient of 
tri-pole mode. Here we plot the ratio of the amplitude of 0(13 and in log 
scale. The dash curves are numerical results, and the solid curves are 
analytical results. The black, red and blue curves are for the case of £1 = 
10"^, £1 = lO"* and £1 = 10"^, respectively. The permeability of ENZ in all 
cases is set as 0.6. The radius of the cylindrical defect is 0.01 m, and its 
dielectric constant is £21 = 16 and its relative permeability is /i2i = 1. In the 
simulations, we set a = h = 0.021 m for the ENZ area. 
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Figure 7 | The magnetic field distribution of ENZ waveguide system with 
a defect for (a) Hi = i and (b) Hi = 1. The working frequency is 
7.615 GHz. 



H2 = zHiY, 



Joik 



/o(M;; 



(28) 



which is consistent with the result found in Ref. 20. In this situation, 
only the monopole mode is excited inside each defect. 

Discussions 

We find some interesting phenomena in ZIM waveguide system 
embedded with defects when total transmission occurs. For 
MIZIM case, additional dipole modes, besides monopole modes, 
could be excited in the defect. For ENZ case, additional higher order 
modes can also be excited if two particular conditions are satisfied at 
the same time. We have shown the underlying physics of why higher 
modes appear. We hope that these phenomena could be observed in 
experiments in the coming future, considering the current experi- 
mental progress on ZIM"' ''''". 



7.61 GHz, as already shown in Fig. 6). Hence, the simulation results 
prove our finding, and the magnetic field is a summation of mono- 
pole mode and tri-pole mode, which should be written as, 

H2 = zHi f^")^^ +ai3/3(fc2iri)cos(30)) . (26) 

\yo(,'c2iJ<ij / 



In order to further demonstrate Eq. (25), we set Is^ii^i) — 0. 
where fc2ii^i is the second root of the fifth order of Bessel function. 
The radius ofthe defect is 0.01 m, its dielectric constant is £21 = 14.44 
and its relative permeability is ^21 = 1. In order to satisfy ]4}^2\^\) ~ 0- 
the working frequency is about 15.5 GHz. To meet Eq. (24), we set 
the effective permittivity and permeability of the ENZ as (,y = 0.001 
and fii = 0.6427 respectively and set a = h = 0.021 m for the ENZ 
area. Fig. 8 is the magnetic field distribution for the new system when 
total transmission happens. The penta-pole mode comes up inside 
the defect when Hi = i, as shown in Fig. 8(a). By changing Hi into 1, 
the monopole mode shows up inside the defect, as shown in Fig. 8(b). 
Besides, we numerically find that the coefficient of penta-pole mode is 
= 1.83;, /J15 = 0 and the coefficient of the monopole mode is 
1/70(^21^1) = 8.46. Therefore, the magnetic field is a summation of 
monopole mode and penta-pole mode, which should be written as. 



H2-ZH1 +ai5/5(fc2iri)cos(5^?) 
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In addition, if only Eq. (24) is satisfied, but none of /„(/c2 jRj) is zero, 
the magnetic field inside each defect is written as, 
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Figure 8 | The magnetic field distribution of ENZ waveguide system with 
adefectfor (a) Hi — /and (b) Hi — 1. The working frequency is 15.5 GHz. 
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